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ABSTRACT 

We present the first galaxy-galaxy weak lensing results using early data from the Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS). These results are based on ^22 deg 2 of i' data. From 
this data, we estimate the average velocity dispersion for an L* galaxy at a redshift of 0.3 to be 137 
± 11 km s _1 , with a virial mass, M200, of 1.1 ± 0.2 xlO 12 h Mq and a rest frame mass-to-light 
ratio of 173±34 hM Q /L Rc Q. We also investigate various possible sources of systematic error in detail. 
Additionally, we separate our lens sample into two sub-samples, divided by apparent magnitude, thus 
average redshift. From this early data we do not detect significant evolution in galaxy dark matter 
halo mass-to-light ratios from a redshift of 0.45 to 0.27. Finally, we test for non-spherical galaxy dark 
matter halos. Our results favor a dark matter halo with an ellipticity of ~0.3 at the 2 a level when 
averaged over all galaxies. If the sample of foreground lens galaxies is selected to favor ellipticals, the 
mean halo ellipticity and significance of this result increase. 

Subject headings: gravitational lensing, dark matter, mass-to-light ratios, galaxy halos, halo shapes 



1. INTRODUCTION 

It is widely accepted that galaxies live in massive 
dark matter halos, but the properties of these halos 
are not particularly well understood. Dark matter ha- 
los around individual galaxies introduce small coherent 
distortions to the shapes of background galaxies. This 
signal, weak gravitational lensing, can be used to in- 
fer properties of foreground galax y dark matter halos 
such as their sizes and shapes (e.g , 
Hudson et al. 1998; Hoekstra et al 



Braincrd et al.11199 
2003; Sheldon et all 



2004t iHoekstra et all 12004 120051: [MandclbauinetalJ 



2006al iHevmans et al.1 l2006bl ). Galaxy-galaxy lensing 
also provides an important link between numerical sim- 
ulations, which model the dark matter very well, and 
other observational techniques, which are restricted to 
studying the (likely) biased luminous matter. The abil- 
ity to connect observed galaxies to the properties of their 
dark matter halos provides important insights into the 
details of galaxy formation. Galaxy-galaxy lensing can 
also be used to constrain alternative theories of gravity. 
If observations such as flat rotation curves are due to a 
modified gravity law then the weak lensing signal at large 
distances from the galaxy center should be isotropic. If 
an anisotropic galaxy-galaxy lensing is observed around 



galaxies then this is a strong evidence for flattened dark 
matt er halos and disfavors an y modified gravity theory 
(e.g. iMortlock fc Turnerll2001t ). 

Other techniques to map the dark matter content of 
galaxies, including extended rotation curves and dynam- 
ical methods, require visible tracer populations and are 
thus capable of probing only the inner regions of galaxy 
halos. Using satellite galaxies to probe the potential 
wells of galaxies offers an interesting alternative (e.g., 
Zaritskv fc Whitdll994l : IPrada et al.ll2003HConrov et al.1 
2005f ). but this technique requires the assumption of dy- 
namical equilibrium and is susceptible to the inclusion of 
interlopers which can bias the results. 

The study of galaxy-galaxy lensing has grown dramati- 
cally since its first detection 10 years ago (Braincrd et al. 
1996), largely due to improved analysis techniques 
and a wealth of wide-field data. The basic re- 
quirement for a weak lensing measurement is wide- 
field imaging of reasonable depth and image quality. 
The most precise galaxy-galaxy lensing measurements 
to date are from t he Red-Sequence Cluster Survey 
(jHoekstra et al.ll2004T) and the Slo an Digital Sky Survey 
(SDSS) (|Mandelbaum et al.ll2006aD which are not partic- 
ularly deep surveys, but cover enormous areas. Measure- 
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merits have been m ade from space in small fields (e.g., 
iHudson et alJll998ft and from large ground-based obser- 
vations (e.g.. ISheldon et alj 12004 I Hoekstra et "all 12004 
iMandelbaum et al.ll2006af) . 

In this study we examine the weak lensing signal from 
a ~22 square degree area of sky from 2 fields of the 
WIDE component of the Canada France Hawaii Tele- 
scope Legacy Survey (CFHTLS) 1 . These data represent 
a small fraction of the total CFHTLS- WIDE area which 
will be covered by the end of the survey in 2008 (approx- 
imately 170 square degrees). This early analysis is based 
on i' data only and is therefore lacking colors and redshift 
information for both the lenses and the sources. Com- 
pared to pre vious large surveys u s ed for galaxy-galaxy 
lensing (e.g.. IHoekstra et al.|[200l ISheldon et all 12004 
the CFHTLS is deeper and therefore we have a higher 
source density and can probe galaxy halos at systemati- 
cally higher redshifts. 

In this paper we will introduce the CFHTLS data (§2) 
and explain the redshift distribution determination and 
the galaxy shape measurements. In §3 we will show the 
first CFHTLS galaxy-galaxy lensing measurements in- 
cluding the determination of the average halo velocity 
dispersion as well as estimates for the average galaxy 
mass and mass-to- light (M/L) ratio. We also investigate 
a number of possible sources of systematic error. We will 
conclude this section by dividing the lens sample into a 
low and high redshift sample in order to look for evolu- 
tion in the galaxy halos. In §4 we will present a measure- 
ment of the halo shapes using anisotropic weak lensing. 
In §5 we will summarize and discuss these results. 

Throughout the paper a cosmology with fi m = 0.3 and 
Qa = 0.7 is assumed. Results are presented in units of 
the hubble parameter, h (the Hubble constant rescaled 
in units of 100 km s _1 Mpc _1 ), which is assumed to be 
1. 

2. DATA 

Canada and France have united to use a large frac- 
tion of their telescope time at the Canada- France-Hawaii- 
Telescope (CFHT) in order to complete a 5-year photo- 
metric survey. The survey makes use of the new Mega- 
Cam instrument at CFHT. The camera is comprised of 
36 separate CCD chips and produces distortion corrected 
1 sq. degree images with superb image quality. The 
survey is divided into DEEP, WIDE and VERY WIDE 
components. The wide component is primarily designed 
for weak lensing studies. These data are well suited to 
the study of galaxy dark matter halos, as will be dis- 
cussed in this paper, but also the study of lensing by 
the large s cale mass distributio n in the universe, cos- 
mic shear (Hockstr a et al.l 12006). Cosmic shear studies 
complement the Type la supernovae analy sis from the 
CFHTLS-DEEP survey ((Astier et alJ 12006( 1 since they 
can both be used to determine cosmological parameters, 
in particular the dark energy equation of state. 

This galaxy-galaxy lensing project makes use of early 
data from the CFHTLS- WIDE survey. The wide data 
will eventually cover ~170 sq. degrees in 5 filters 
(it* ,g' ,r' ,i' ,z'). The observations are divided into 4 large 
patches which are well-separated in right ascension. Each 
patch is located far from the Galactic Plane in order to 

1 http:/ /www. cfht.hawaii.edu/Science/CFHLS/ 
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Fig. 1— CFHTLS observations of Wl and W3 fields. The 
effective area is approximately 22 deg 2 from 3f pointings (f9 deg 2 
in Wf and f 2 deg 2 in W3 ). Bright stars and bad pixels are masked. 

minimize extinction and contamination from bright stars. 

The early data used in this paper study is based on 
31 pointings taken in the i'-band and covers an effective 
area of roughly 22 deg 2 . The data are from the Wl and 
W3 fields which will total 72 and 49 deg 2 respectively 
at the end of the survey. Wl is centered at 2 h 18 m 00 s 
-7°00'00" and overlaps with the XMM-LSS field. W3 
is centered at 14 h 17 m 54 s +54°30'31" and overlaps with 
the Groth strip field. Figure Q] shows the coverage of 
the data used in this paper, totaling 19 deg 2 from Wl 
and 12 deg 2 from W3. The images are masked to avoid 
bright stars, diffraction spikes and bad pixels. The source 
density is ~20 galaxies per sq. arcminute. These data 
were obtained in the 2003B, 2004A and 2004B observing 
semesters. The median seeing for the images used in this 
project is 0.76", however, more recent CFHTLS obser- 
vations have improved image quality. The image quality 
degrades towards the edge of the images but only data 
with sub-arcsecond seeing are included in this analysis. 
More de tails on the data use d in this analysis can be 
found in IHoekstra et all (|2006( 1. 

2.1. Images and Catalogs 

The i' images used in this analysis are provided to 
the Canadian and French Astronomical communities via 
the Canadian Astronomical Data Centre (CADC). The 
basic data reduction pipeline, Elixir, is used to provide 
flat-fielded, de-biased images as well as photometric zero- 
points and an initial astrometric solution. The CFHTLS 
data are complicated by the fact that there are large 
dithers between some of the images. Therefore, some ob- 
jects may appear on different chips in the individual ex- 
posures. This can lead to complicated PSF anisotropies, 
which are difficult to model and correct for. With this 
early analysis all such "multi-chip" data are excluded. 
Each chip is stacked separately with data only from that 
chip. This comes at a cost of losing approximately 20% 
of the area but ensures accurate PSF anisotropy correc- 
tion. Using stacks where all of the data comes from 
only one chip is a conservative approach and since the 
Megacam PSF is smooth from chip to chip future anal- 
yses will not impose this restriction. As described in 
IHoekstra et"all (2006) any images with poor seeing (> 1 
arcsecond) are excluded from the stacks. Additionally, 
the images are inspected by eye to mask out areas con- 
taminated by bleeding stars, diffraction spikes or other 
cosmetic defects. The stacks for this CFHTLS data were 
created using the SWarp routine 2 . The stacks are created 
from 6 or 7 individual images of 620 second exposures, 

2 http:/ /terapix.iap.fr/rubrique.php?id_rubrique=49 
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thus each stack is >1 hour and reaches a depth of ~25 
in i' . 

A critical aspect of weak lensing measurements is accu- 
rately determining the shapes of faint objects. The abil- 
ity to do this with stacked images requires very accurate 
astrometric solutions so that no spurious shape distor- 
tions are introduced in the stacking procedure. For this 
reason it is necessary to improve upon the basic astrom- 
etry from the Elixir processed data and include higher 
order corrections. The images and catalogs used in this 
galaxy-galaxy lensing analysis are the same as those used 
in the first cosmic she ar analysis of the CFHTLS- WIDE 
(jHoekstra et al.ll2~006| ) . The catalog details can be found 
in the shear paper, but we will outline the major points 
briefly here. 

The basic technique for correcting the image astrome- 
try is as follows: 

• Retrieve a red image from the second generation 
Digital Sky Survey (POS II) for each pointing 



• Use SExtractor (jBertin fe Arnoutslll996t ) to gener- 
ate a high density catalogue of stars 

• Match astrometric catalogue to each MegaCam im- 
age, averaging the positions of objects which ap- 
pear on more than one chip (due to the large dither 
offsets) 

• Use this master catalog to generate an astrometric 
solution for each chip 

Following the above procedure a stack can be created 
using all of the input images with the new astromet- 
ric solution applied. As was discussed above the stacks 
were created for each chip individually, and include only 
data from that chip, so that the PSF pattern could be 
corrected properly. This reduces the total area but sim- 
plifies the PSF corrections. 

2.2. Galaxy Shapes 

The carefully stacked i' images were processed wit h 
the peak finding algorithm of KSB (Kai ser et al.l fl995). 
Objects which were detected to be more than 5cr above 
the sky background were added to the source catalogue. 
The catalogue is then cleaned to include only objects 
larger than the PSF (points to the right of the stellar 
locus, identi fied in a plot of magni tude versus half-light 
radius as in [Hoekstra et all (|2002f )). These objects are 
then analyzed in more detail in order to determine their 
apparent i' magnitudes and half-light radii, as well as 
their shape parameters. The shape is defined by two 
polarization "vectors" (equations [T] and [2]) : 
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(2) 



where Qij are weighted with a Gaussian function that 
scales with galaxy size 



2.2.1. Shape Corrections 

The measured shapes must be corrected for dis- 
tortions such as the effects of seeing and PSF 
anisotr opy. This is done following the tech niques out- 
lined inlKaiser. Squires, fc Broadhurstl (fl995f) (KSB) and 
iHoekstra et al.l (|1998|) . A sample of stars is found in the 
images and they are used to characterize the seeing and 
any anisotropy in the PSF. The shapes of the stars are 
fit with a second order polynomial for each chip in the 
CFHTLS stacks individually (that is to say, the stars are 
fit in the final stacked image, but each stacked chip is 
handled individually). The shapes of the stars are then 
used to correct the shapes of all of the galaxies. The 
catalogues include the position, magnitude, shape, error 
and P 7 information, where P 7 is the pre-seeing shear 
polarizability and can be calculated fro m the images of 
stars and galaxies (Hoekstr a et al.|[l998ft . 

The procedure outlined here to correct the galaxy 
shapes was recently tes ted in the Shear TEsting Program 
(iHevmans et a l. 2006a), and with more realistic PSFs in 
iMassev et al.1 ( 20071 ) . The results indicated that the pro- 
cedure followed to create the catalogues in this analysis 
can reliably be used to measure weak shear to an accu- 
racy within ~2 %. It is impor t ant to not e that the mor- 
hologi es in the lMassev et al.1 (|2007h and lHeymans et all 
2006a) papers are simplified and do not consider possi- 
bilities such as shear and noise that vary a s a function of 
positio n. The simulated images used in the lMassev et al.l 
(2007) analysis are based on the best available imaging 
from Suprime-Cam and thus do not have the identical 
signal-to-noise or PSF size as the CFHTLS images, there- 
fore the quoted 2% shear accuracy should be taken as a 
rough estimate. However the lack of B-m odes in the cos- 
mic s hear analysis using these catalogues (jHoekstra et al.l 
2006) is a powerful indicator that the level of system- 
atics arising from PSF anisotropy is very small. Fur- 
thermore, galaxy-galaxy lensing is not very sensitive to 
PSF anisotropy as we average over many pairs of galaxies 
which have random directions of PSF anisotropy. 

2.3. Redshift Distribution 

The final CFHTLS- WIDE data set will include pho- 
tometric redshifts for every faint source, but at present 
we have only single filter data available. For the analysis 
presented here, we select a sample of lenses and sources 
on the basis of their apparent i' magnitude. We define 
galaxies with 19<z'<22as lenses, and galaxies with 
22.5 < i' < 24.5 as sources which are used to measure 
the lensing signal. The faint end of the source distribu- 
tion is chosen to ensure high signal-to-noise shape mea- 
surements and corresponds to the peak in the magnitude 
distribution (Figure^. This selection yields a sample of 
^2xl0 5 lenses and ~1.3xl0 6 sources. These catalogues 
are used to generate 3.7xl0 7 lens-source pairs within a 
projected radius of 2 arcminutes of the stacked lenses. All 
lens-source pairs within 7 arcseconds of the host galaxy 
are eliminated from the catalogue since their shape mea- 
surements are likely compromised by light from the host 
lens. The angular scale of 7 arcseconds to 2 arcminutes 
corresponds to a physical scale of ~25-500 h _1 kpc at the 
median redshift of the lenses. Figure [2] shows the magni- 
tude distribution for the entire sample of galaxies in the 
Wl and W3 fields used. 
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Fig. 2. — Lens and source magnitude distributions. This his- 
togram shows the distribution of i' magnitudes for all galaxies in 
the Wl and W3 fields used in this analysis. The hatched region 
divides the lenses and the sources, and the vertical line at i'=24.5 
shows the upper magnitude for the sources used. 

The lensing signal for an isothermal sphere is a func- 
tion of (f3), the average ratio of the angular diameter 
distances between the lens and source, D^s, and between 
the observer and the source, D5 as follows 



0. 



D 



LS 
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Therefore, in order to interpret the detected shear mea- 
surements it is necessary to know the redshifts of both 
the lenses and the sources. If the precise redshifts are 
not known for each object then at least their distribu- 
tions must be understood in order to convert shear mea- 
surements into physical properties such as velocity dis- 
persions and halo masses. The CFHTLS data used here 
were taken in a single band and so the redshift distribu- 
tions of the sources and lenses must be estimated. The 
shear can only be estimated in projected angular bins 
and not physical units such as kiloparsecs. Therefore, 
the lensing signal for a distant galaxy is measured on a 
much larger physical scale than for the case of a closer 
galaxy. The mixing of scales complicates the interpre- 
tation of the results, but we can still learn about the 
average properties of halos. 

The lens catalogues contain relatively bright galaxies, 
and therefore the distribution of the redshifts is quite well 
understood from pr evious studies, s uch as the CNOC2 
field galaxy survey (|Yee et al.l [2000). It is much more 
difficult to estimate the redshifts of faint sources, since 
they are generally too faint for spectroscopic redshift de- 
termination. The redshifts of faint galaxies are generally 
estimated using photometric redshifts. A spectroscopic 
study of the Hubble Deep Field (HDF) to depths of -24 
in R c showed that the spectrosc opic redshifts agree d well 
with the photometric redshifts ([Cohen et al.ll200"o[ ). 

The source and lens redshift distributions used 
in this analysis can be seen in Figure [3J The me- 
dian lens redshift is ~0.4 and the median source 
redshift is —0.9. The N(z) distribution for the 
lenses was estimated using the functional form of 
iBrown. T aylor. Bacon . Gray. Dye. M ciscnhcimcr. & Wolf 



Fig. 3. — N(z) distribution. The lens and source distributions 
used in this analysis are shown with the solid lines. The sources 
have a median redshift of ~ 0.9 and the lenses have a median 
redshift of ~ 0.4. The dashed line shows the redshift distribution 
of the sources based on the photometric redshifts from the HDFs 
HFernan dcz-Soto et al. 1999). The sources are at lower redshift 
using the HDF photometric redshifts, but the value of (/3) is only 
~5% smaller. 

(pOOl based on COMBO-17 data. The COMBO-17 
survey uses a combination of 17 wide and narrow band 
filters to measure very accurate photometric redshifts in 
a 1 sq. degree patch. 
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where z* is related to median redshift of the distribution 
by z* = z m /1.412, and z m is a function of the apparent 
magnitude. 

The source redshift distribution used in the anal- 
ysis was estimated fro m photometric red shifts in the 
CFHTLS-DEEP fields filbert et al.ll200l . The func- 
tional form of the N(z) sources was as follows 



N(z) s 



where 
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The values for A, lu and zq were calculated using the 
magnitude c uts of our source catalogues applied to the 
lllbert et al.l (|2006l ) data. 

The N(z) distributions can then be used to estimate 
the angular diameter distances to the lenses and sources. 
Assuming a standard ACDM cosmology this results in 
a (f3) of 0.49. If we alternatively used a source dis- 
tributi on based on the Hubble De ep Fields North and 
South (|Fernandez-Soto et al.l Il999h then we find a ((3) 
of 0.46±0.02, where the error is from the field-to-field 
variations in the HDFs, and the finite number of galax- 
ies. The P estimates are weighted in the same manner 
as the shear estimates as will be discussed in the next 
section. The difference is source redshift distributions 
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between the CFHTLS-DEEP fields and the HDFs is ob- 
vious in Figure [3l The sources are at slightly higher 
rcdshifts when using the CFHTLS photometric redshifts 
than with distribution estimated from the HDF. The vis- 
ible discrepancy between the photometric redshift distri- 
butions from the CFHTLS and the distributions based 
on the HDF is not statistically significant when one ac- 
counts for uncertainties in photometric redshift distribu- 
tions d ue to cosmic varianc e affec ting small fields (see 
Fig 1 of Ivan Waerbeke et alj (|2006f ) for a dramatic illus- 
tration of this in the HDF). The different redshift distri- 
butions only change the (0) value by ~5%. 

3. ANALYSIS 

The tangential shear signal must be fit with an as- 
sumed halo mass model in order to extract physical prop- 
erties of the halo such as velocity dispersion and mass. 
A common mass model assumed is the isothermal sphere 
for which the tangential shear is proportional to the Ein- 
stein radius and hence to the velocity dispersion squared 
(equations [7] and [SJ . 
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In this study the lenses are stacked together and the 
sources that lie within a projected radius of 2 arcminutcs 
are divided into angular bins. The component of their 
shape tangential to the lens center is determined and av- 
eraged in each bin. Each shear calculation is weighted by 
the error in the shape measurement as described in Hoek- 
stra et al. (2000). Galaxy-galaxy lensing measurements 
also have a convenient built-in systematic test. If the 
tangential lensing signal detected is due to gravity then 
it should vanish if the source images are rotated by 45 
degrees. The tangential and "cross-shear" for the entire 
sample are plotted in Figure 2] The best fit isothermal 
sphere has an Einstein radius of 0'.'247±0'.'020. A lensing 
signal is detected at high significance (> 12cr). The cross- 
shear measurement is consistent with 0, as expected, and 
therefore the tangential shear signal is interpreted as be- 
ing caused by weak lensing from an isothermal sphere 
potential. 

The tangential shear measurements can also be fit with 
other dark matter profiles, such as the Navarro , Frenk 
and White (NFW) profile l|Navarro et al.lll996l ). given 

by 



p(r) 



ScPc 



(r/r«)(l + r/r s 



(9) 



where p c is the critical density for closure of the universe. 
The scale radius, r s , is defined as r2oo /cnfw where 
cnfw is the dimensionless concentration parameter, and 
5 C is the characteristic over-density of the halo. The tan- 
gential shear equat ions for a NFW halo can be fo und in 
iBartelmannl (|2006l ) and lWright fc Brainerdl (|200Cf >. The 
NFW density profile has 2 free parameters, but in this 
analysis we will consider the concentration parameter 
to be fixed, obeying the relation between virial mass 
and co ncentration found in simulations by Bullock et al. 
(|2001f h The best fit NFW profile, corresponding to a 
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Fig. 4. — (a) The ensemble averaged tangential shear as a 
function of radius around a sample of CFHTLS galaxies with 
19 < i' < 22. The best fit isothermal sphere, shown with the solid 
line, yields an Einstein radius of 0"210±0'.'017, corresponding to a 
velocity dispersion of 132±10 km s — 1 . The dashed line represents 
the best fit NFW profile, corresponding to a halo with an r2oo of 
150 h — 1 kpc. (b) The signal when the source images are rotated 
by 45° . Note the difference in scales between the top and bottom 
figures. No signal is present as expected if the signal in (a) is due 
to gravitational lensing. 

halo with r 200 of 0.15h _1 Mpc (M 200 of 7.6x lO 11 !^^©) 
is shown as the dashed line in Figure [4] 

3.1. Systematics 

In order to verify the validity of our shear signal we 
have performed a number of systematic tests. At the 
outset we are confident that the catalogue used in this 
analysis is free of major systematics because there were 
no significant B-modes found in the cosmic shear analy- 
sis using the same catalogue (|Hoekstra et al.l l2006). The 
initial systematic test is to calculate the cross-shear, as 
discussed above. We have also considered the fact that 
some lenses and sources are physically associated. We 
can estimate the level of this association by calculating 
the number density of galaxies in each bin, which de- 
creases with radius. We have estimated the fractional 
excess of background sources around lenses, fo g , and in- 
creased the shear signal by 1 + fo g to compensate. The 
background density of sources was estimated by count- 
ing the number of sources in bins around random lens 
positions. The random lens catalogue was the same size 
as the lens sample used in this analysis. Our best fit 
fog = 0.82r~ - 63 , where r is the radius in arcseconds, 
is similar t o what has been fou nd in previous measure- 
ments (e.g. Hoekstra et al. 2004) The fractional excess is 
a strong function of radius and thus only the inner an- 
gular bins are influenced by physically associated lens- 
source pairs. This correction has been applied to all 
quoted values in this paper (e.g., 6e, cr), and was cal- 
culated for each lens sample separately. 

A further test of the signal can be made by measur- 
ing the tangential and cross-shear around random points 
in the field. The results from performing this test are 
demonstrated in Figure [5] and are consistent with no de- 
tected signal, as expected. The same number of random 
points were used as the number of foreground lens galax- 
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Fig. 5. — The ensemble averaged tangential (filled circles) and 
cross (open squares) shear as a function of radius around a sam- 
ple of random points in the CFHTLS fields. Both the signals are 
consistent with zero as expected. Note that the cross-shear points 
have been slightly shifted so that the error bars are clearly visible 

ies. 

It is also important to examine the influence of intrinsic 
alignments. Without redshift information for our lenses 
and sources it is possible that physically associated lens- 
source pairs arc included in the analysis. If these pairs 
are aligned then the tangential shear signal will be con- 
taminated. Recent evidence has suggested that satellite 
galaxies' maj or axes preferentially point towards their 
host galaxies (jAgustsson fc Brainerdl[2006h which would 
decrease the galaxy-galaxy lensing signal. If we adopt the 
most pessimistic view and consider a very large contami- 
nation in the inner 250 kpc due to intrinsic alignment, we 
can compensate by boosting the tangential shear signal. 
If we apply a 40% boost to the inner 3 bins in Figure 
[4] (corresponding to roughly 250 kpc at the redshift of 
our lenses), adopting the worst case from Agustsson & 
Brainerd (25-40% suppression of shear on scales less than 
250 kpc), then the best fit isothermal sphere yields an 
Einstein radius which is 7% larger than we have quoted 
above. This small difference is due to the fact that the 
outer 3 data points, with smaller errorbars, constrain 
the isothermal sphere fit much more than the 3 inner 
bins. This is an upper limit on the contamination from 
intrinsically-aligned satellites. 

Alternatively, we can also estimate the influence of in- 
trinsic alignment by repeating the analysis after dividing 
the distant sources into 2 sub-samples based on their ap- 
parent magnitudes, while keeping the lens sample iden- 
tical. After re-scaling by the new j3 values the shear 
signal should be the same since the lens sample has not 
changed. If, however, there is significant intrinsic align- 
ment, then the signal from the brighter sources (lower 
redshift) will be suppressed. We have carried out this 
test and found that the velocity dispersion derived from 
the shear signal computed using the brighter sources is 
indeed suppressed by <10%. This suppression due to 
intrinsic alignments for the entire sample would there- 
fore be <10%. We have not applied any correction for 
intrinsically-aligned sources, instead we include this as 



part of our total systematic error budget. 

One final potential systematic that should be discussed 
is possible contamination from stellar objects in the fi- 
nal source catalogue. Moderately bright stars are easily 
identified as they lie in a tight vertical locus in plot of 
magnitude versus size. Bright stars saturate and their 
measured sizes increase, but they can also be easily re- 
moved from the catalogues. There is some overlap be- 
tween faint stars and galaxies which can dilute the shear 
signal, but by extrapolating the number counts of stars 
to faint magnitudes it is clear that stars will contribute 
at most a few percent to the total number counts. In 
addition, we consider only sources larger than the stellar 
PSF which cleanly sep arates stars and galaxies to i 1 ~ 24 
(|Hoekstra et al.|[2002T ) . Furthermore, in this analysis we 
down- weight faint sources which significantly reduces the 
contamination from stars (Heymans et al. 2006). The 
weight for sources fainter than %' of 24 is very low, thus 
the level of stellar contamination is estimated to be <1%. 
It is important to note that for the discussion of halo 
shapes below that this potential source of systematic er- 
ror is irrelevant because stars are randomly distributed 
and the measurement is a relative one. 

Based on the above considerations we estimate that 
the total systematic error in our measurements is at most 
10% (1-cr level). The exact level of possible systematic 
errors is difficult to estimate as the error distributions 
for the various systematics are not known. Throughout 
the paper the quoted errors are always statistical errors 
derived based on the errors in the shape measurements 
and the number of sources in each bin. 

3.2. Velocity Dispersion 

The best fit isothermal sphere yields a velocity disper- 
sion of 132±10 km s" 1 (equations [7] and [8|) . This mea- 
surement takes into account the fact that some lenses 
and sources will be physically associated. We can esti- 
mate the level of this association by calculating the num- 
ber density of galaxies in each bin which decreases with 
radius. 

The velocity dispersion depends on the sample of lenses 
used and must be scaled in order to compare to other 
results. This can easily be done by assuming a scaling 
relation between luminosity and velocity dispersion as 
shown in equation 1101 



where o* is the velocity dispersion of an galaxy. The 
scaling factor a is generally assumed to be 3 or 4, mo- 
tivate d by the observations of the Faber- Jackson rela- 
tion (|Faber fc Jacksonl Il976f) . for example. We based 
our L » galaxy on th e CNOC2 luminosity function re- 
sults (|Lin et al.lll999f ). Lin et. al measure L* at z= 0.3, 
which after correction to the ACDM cosmology assumed 
here gives 1.3 x 10 10 fe~ 2 L^H . We u sed the color and k- 
corrections from lFrei fc Gunr] (|1994f ) to convert our mag- 
nitudes to R c and evaluate the luminosity of our lenses. 
All luminosities are calculated assuming the proportion 
of early and late types of galaxies found in the CNOC2 
redshift survey, since the redshifts and brightnesses of 
the two samples are comparable (28% early types, 24% 
intermediate types, 47% late types). It is also important 
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to note that the average luminosity is calculated using 
the same weights used in the shear analysis. 

The Einstein radius for an L* galaxy is given by 



4tt oj 

C 2 L 2Ja 



(PL 



2/a\ 



(11) 



where <7* is the velocity dispersion for an L* galaxy and 
6e is in radians. We estimate the average luminosity 
for our lens galaxies to be (L) = l.lxl0 10 /i~ 2 L Q fl c . The 
scaled velocity dispersion can now be estimated for differ- 
ent assumed a values. The results of scaling the observed 
velocity dispersion for our sample of lenses to a typical 
L* galaxy are summarized in Table [1] 

Note that, in principle, the tangential shear signal is 
due not only to the lens galaxy in question (the 1-halo 
term, assuming the galaxies are central galaxies), but 
also other mass clustered with the primary galaxy (the 
2-halo term). In practice, though the velocity disper- 
sion is estimated from the best-fit isothermal sphere out 
to a radius of ~130", which at these redshifts corre- 
sponds to a physical scale of much less than 1 h _1 Mpc 
which should still be dominated by the 1-halo term 
(|Zehavi et alj |2004) . This is not strictly true for satellite 
galaxies where the contribution to the lensing s ignal from 
the h ost halo would kick-in on smaller scales (|Yoo et al.l 
2006). However, as the shear profile is well fit by both 
Singular Isothermal Sphere and NFW profiles this is not 
a significant concern. 

A galaxy-galaxy lensing analysis of the COMBO-17 
data by Kleinheinrich et al. (2005) estimated the value 
for <7* for various scenarios, including when the redshift 
is known for every lens and every source, and when there 
is no redshift information available. The latter case can 
be compared with our results. They obtained a mean 
er» of 156^3Q km s _1 . Hoekstra et al. (2005) used data 
from the Red-Sequence Cluster Survey (RCS) to estimate 
the properties of galaxy dark matter halos using galaxy- 
galaxy lensing without redshift information. They found 
a value for er* of 136±8 km s — 1 . It is important to note 
that the estimate for L» is slightly different in each of 
these studies, but the results agree within the errors. 

3.3. Halo Masses 

The total extent and mass of dark matter halos can be 
estimated ass uming a mass model for the ga l axy ha los as 
suggested by iBrainerd. Blandford. fe Smaill (|1996h with 
a density profile 



p(r) 



2irGr 2 (r 2 + s 2 ) 
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where s is a measure of the truncation scale of the halo. 
This profile is an isothermal sphere at small radii with a 
cut-off at large radii, characteri zed by the scale s, whi ch 
scales with velocity dispersion (Schneider & Rix 1 997( 1 
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Assuming this truncated isothermal sphere halo model, 
the mass enclosed within a sphere of radius r is 
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which, because of the truncation, results in a finite mass 
(jHoekstra et al.ll2004ft 
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The truncation scale can be assumed to be constant for 
all halos, in which case the M/L ratio would be oc L 1 / 2 , 
if L oc a 4 . This is the favoured model in recent analyses 
fe.g. lGuzik fc Selialdl200H IHoekstra et al.ll2005D . Alter- 
natively, it can be assumed that M /L is constant for all 
galaxy halos in whi ch case s oc a (jBrainerd et al.| [l996; 
iHudson et al.l lT998). We do not explore the various pos- 
sible M/L scalings in this paper as we assume a value 
for the truncation radius from the literature and do not 
directly fit a truncated isothermal sphere model. 

Assuming the trunc at ion radius found by 
iHoekstra. Yee. fc Gladdersl (j2004l ) for an L* galaxy 
of 185±30 kpc we estimate the total mass of our L* 
galaxy to be 2.2 ± 0.4 xlO 12 br x M. & if Loc a 4 , and 
2.4±0.5 xlO 12 /j^Mq if Lex a 3 . The results are in 
good agreement with the results from the Red-Sequence 
Cluster Survey (jHoekstra et al.| [2004). which found total 
halo mass to be 2.7±0.6 xlO 12 /i _1 M . We estimate 
the virial mass, M2oo, of our galaxy halos, assuming an 
isothermal sphere model, to be 1.1±0.2 xlO 12 h~ 1 M Q , 
where M200 is the mass inside the virial radius r2oo 
defined as the radius where the mass density is 200 
times the critical density. This is in good agreement 
with the results from RCS where the virial mass of 
L* halos was estimat ed at 11.7 ± 1.7 x 10 11 h~ 1 M Q 
(jHoekstra et al.1 l2005h . Our results are also in good 
agreement with results from the lower redshift SDSS 
sample where the virial ma ss of L* halo was found t o 
be 13.2i^;° x 10 u h- 1 M o (iMandelbaum et all l2006ah . 
In order to make this comparison we chose the central 
virial mass found for <L /L»>=1.1 (see Table 3 of 
IMandelbaum et aTl (|2006a)) and scaled this to L*. We 
then scaled the SDSS from an NFW with a virial radius 
defined using p=180p to a singular isothermal sphere 
profile with a virial radius defined by p—200p cr it- . The 
GEMS team has estimated the virial mass of halos fo r 
a sample at higher redshift (jHevmans et al.l l2006bf) . 
They estimate that the virial mass for an L* galaxy at 
a redshift of - 0.65 is 14.1±|;| x 10 u h- 1 M Q . For this 
comparison we once again transformed from an NFW 
to a singular isothermal sphere profile, for an L* galaxy, 
using our definition of the virial radius. Our measured 
velocity dispersion for an L* galaxy is also in agreement 
with results from fundamental plane measurements 
(jSheth et al.ll2003h . 

We can use the luminosity of an L* galaxy together 
with the mass estimates above in order to calculate a 
typical M/L ratio. This leads to a M/L ratio for an L* 
galaxy of 173±34 hM Q /L Rc 0, assuming a — 4. 

3.4. Evidence for evolution? 

The CFHTLS data are ideally suited to looking for 
evolution in dark matter halos with redshift. The high 
signal-to-noise of our primary result permits the division 
of the lenses into sub-samples at different redshifts. With 
this early data there is sufficient signal-to-noise to divide 
the foreground lens samples into multiple bins, however 
without redshift or color information this proves diffi- 
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Fig. 6. — (a)The tangential shear signal for two samples of fore- 
ground lenses, (b) The cross-shear signal when the source images 
are rotated by 45°. The faint lenses (indicated by open squares) 
have i' > 20.5 while the bright lenses (indicated by filled circles) 
have i' < 20.5. 

cult. Nevertheless, as was indicated earlier we can divide 
the lenses based on their observed i' magnitudes. The 
brighter lenses will be on average at a lower redshift while 
faint lenses will be, on average, at higher redshift. We 
decided to apply a magnitude cut at 20.5 in i'. This leads 
to a "bright lens sample" containing ~4x 10 4 galaxies 
and a "faint lens sample" containing ~ 1.5 x 10 5 galax- 
ies. This selection results in a low redshift sample with a 
median redshift of 0.27 and a high redshift sample with 
median redshift of 0.45. 

We repeated the tangential shear analysis described in 
the previous section and the results are shown in Fig- 
ure [6] There is a striking difference in the tangential 
shear profiles for the two lens sub-samples. However, as 
indicated in equation[7l the velocity dispersion is also de- 
pendent on the parameter (3. The source catalogues used 
are the same as in the analysis of the entire sample, and 
therefore so is the source redshift distribution. However, 
these two samples of lenses have very different redshift 
distributions, thus their average (3 values will also be dif- 
ferent. The faint lens population has a (/3)= 0.44 while 
(/3)=0.67 for the bright lenses. 

The differences in (3 can explain most of the offset be- 
tween the 2 shear profiles in Figure O however it is still 
interesting to scale the velocity dispersions to a typical 
L» galaxy to see if any real difference can be seen in the 
dark matter halos of the two lens samples. The results 
are shown in the Table [TJ We do not see conclusive ev- 
idence for evolution in galaxy dark matter halos in this 
redshift range, but this is not so surprising due to the lack 
of photometric or spectroscopic redshift information to 
clearly divide the lens populations. 

4. HALO SHAPES 

An important insight into the nature of dark matter 
comes from the shapes of dark matter halos. Dynam- 
ical measurements can be used to trace out the halo 
shapes on small scales, but they can not be used on 
larger scales where there are no visible tracers. Nu- 
merical simulations of CDM indicate that dark matter 



halos s hould be flattened, and mo r e often prolate than 
oblate (iDubinski fc Carlberd ll991t ISpringel et all 12004 
lAlkood et alB ooe). Simulations of self-interacting dark 
matter produce more spherical halos (|Dave et al.l 2001). 

Alternative theories of gr avity such as M odified New- 
tonian Dynamics (MOND) iMiigroml (|2002j ) attempt to 
explain astrophysical observations by modifying gravity 
rather than invoking non-baryonic dark matter. Non- 
relativistic theories of modified gravity can not provide 
predictions for lensing measurements, and therefore can 
not be tested by lensing observations. However, there 
is now a candid ate relativis tic modified gravity theory 
as presented by iBekensteTnl (|2004| ) which can be used 
as an alternative to dark matter to explain relativistic 
phenomena such as lensing. Thus far this theory ap- 
pears to match m any observations (jSkordis et al.l (20061 : 
IZhao et al"1l2006[ ). One interesting test for this theory 
is the inferred sh apes of galaxy dark matte r halos from 
weak lensing (e.g. lMortloc k"fc Turner! feOO If ). 

Modified gravity theories predict that the lensing sig- 
nal is due to the observed luminous material, and thus 
any anisotropy in the lensing measurement is due to the 
anisotropy in the distribution of gas and stars. There- 
fore, on small scales one would expect the lensing signal 
to be anisotropic since galaxies themselves are, but on the 
large scales probed by galaxy-galaxy lensing one would 
expect a nearly isotropic signal since there is no luminous 
material present at large radii from the galaxy center 
(provided that the galaxy is isolated and not in a group 
or cluster). If a highly anisotropic signal is detected at 
large radii this provides an interesting constraint on mod- 
ified gravity theories, and provides supporting evidence 
for dark matter theories. The predictions from MOND 
of spherical halos assumes that the lens galaxies are iso- 
lated, and not in groups or clusters where there could be 
light contamination on large scales. 

Galaxy-galaxy lensing measurements have gener- 
ally assumed that dark matter halos are spher- 
ical, but a recent gala x y-gala xy measurement by 
iHoekstra. Yee. fc Gladdersl (|2004h detected a signifi- 
cant flattening of dark matter halos. This result was 
not observed in the latest analysis of SDSS data by 
iMandelbaum et ail (|2006b). However, Mandelbaum et 
al. did find marginal evidence for flattened halos when 
they restricted their lens sample to red galaxies. 

One important note is that all such measurements rely 
on the assumption that the mass distribution of the halo 
is aligned with the light distribution of the galaxy. If 
the halo flattening is not correlated with the light profile 
orientation then this simple measurement is much more 
difficult to interpret, and the flattening signal is likely 
systematically suppres sed. 

Brainerd & Wright (2000) discussed how to measure 
halo shapes by looking for an anisotropic galaxy-galaxy 
lensing signal. One simple approach to try and detect 
if galaxy dark matter halos are non-spherical is to test 
whether the tangential shear signal is different along the 
semi-major and semi-minor axes of the visible galaxy. 

The analysis described in the previous sections was 
repeated for the galaxy lenses, this time dividing the 
sources into those within 45 degrees of the semi-major 
axis and those within 45 degrees of the semi-minor axis 
(see the schematic in Figure [7]) . The tangential shear 
results can be seen in Figure [U The signal from the 
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TABLE 1 

Properties for the various lens samples. The quoted errors do not include mass model uncertainties. 
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Fig. 7. — Schematic of anisotropic shear. If galaxy halos are not 
spherical then there should be a difference in the tangential shear 
signal coming from the regions labeled with an A versus those 
labeled with a B. 

two angular bins of sources is very similar. The best fit 
isothermal spheres yield Einstein radii of 0'.'211±0'.'021 
and 0'.'270±0'.'031, corresponding to velocity dispersions 
of 122±12 km s _1 and 138±15 km s _1 , respectively. 

We calculated the ratio of (7} minor to (7) major out to 
70" (which corresponds to ^250 h _1 kpc, for comparison 
to Brainerd & Wright) and the results are plotted in 
Figure[9l The best fit shear ratio is 0.76±0.10, indicating 
a ^2<7 detection of non-sphericity for dark matter halos. 
Examination of the minor-to-major axis shear ratio as 
a function of radius in Figure 1 of Brainerd & Wright 
(2000) suggests that a ratio of 0.76 corresponds to a halo 
cllipticity of ~0.3. 

It is important to measure the anisotropic weak lens- 
ing signal for a sample of well-understood lenses. A way 
to improve the measurement of anisotropic weak lensing 
is to select a sample of lenses with a well defined semi- 
major axis direction. Therefore it is perhaps prudent to 
choose only the brighter lenses, or discard those galax- 
ies with noisy shape measurements, or with very little 
ellipticity. If a lens has an extremely small measured el- 
lipticity, then determining an accurate semi-major axis 
location is very difficult and prone to errors. We choose 
to repeat the anisotropic weak lensing signal for all lenses 
with an ellipticity >0.15 which corresponds to an axis ra- 
tio <0.8. The results are shown in Figure ITOl The best 
fit shear ratio is 0.56±0.13 which according to Brainerd 
&: Wright favors halo model with an cllipticity of ~0.5, 
although directly comparing to Brainerd & Wright is not 
strictly correct as their results are based on simulations 
including the full distribution of galaxy shapes. 

In addition, it is of interest to divide the lens popu- 




Fig. 8. — Mean tangential shear for sources close to the semi- 
major axes (filled circles) and close to the semi-minor axis (open 
squares). The best fit isothermal sphere for the sources within 45° 
of the semi-major axes, as indicated by the dashed line, yields an 
Einstein radius of 0'.' 228±0'.' 026 corresponding a velocity dispersion 
of 138±15 km s _1 . The best fit isothermal sphere for the sources 
within 45° of the semi-minor axes, as indicated by the solid line, 
yields an Einstein radius of 0'.'184±0'.'018 corresponding a velocity 
dispersion of 122±12 km s . The cross-shear is consistent with 
as expected. 
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Fig. 9. — Ratio of mean shear experienced by sources closest to 
minor axes of a foreground lens to those experienced by sources 
closest to the major axes. The weighted average shear ratio is 0.76 
±0.10 favoring a halo with an ellipticity of 0.3. 
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Fig. 10. — Ratio of mean shear experienced by sources clos- 
est to minor axis of lenses to those experienced by sources closest 
to the major axes as in Figure [9] only with the roundest lenses 
(c< 0.15 excluded). Lenses which appear circular on the sky have 
more poorly determine semi-major axes positions and therefore 
have been excluded from this plot. The weighted average shear 
ratio is 0.56±0.13. 

lation into red and blue (or early- and late-type) sub- 
samples and to focus mainly on the early type galaxies, 
since they will not be so susceptible to inclination effects 
(|Novak et al. 2006). Without color or morphology infor- 
mation this was a difficult task with our data, but by 
selecting a sample of galaxies with an axis ratio, b/a, 
between 0.5 and 0.8 we are se lecting mostly early type 
galaxies (jAlam fc Rvdenll2002| ). The results from apply- 
ing this cut in axis ratio to our lenses is shown in Figure 
111! Using preferentially early type galaxies yields a shear 
ratio of 0.61±0.14. The favors a halo with ellipticity of 
~0.4 and represent a more significant detection of non- 
spherical halos than when using the entire lens sample. 
It is important to note that in this analysis we make the 
assumption that the galaxies we are using are isolated, 
so a comparison to MOND is perhaps not completely 
justified. However, the lensing signal in this analysis is 
measured on scales dominated by the contribution to the 
" 1-halo" term from the matter distribution around the 
galaxy itself rather than from the host halo of any satel- 
lite galaxies within groups or clusters. 

A recent analysis of isolated galaxies in the SDSS in- 
dicated that the number density of satellite galaxies is 
higher near the major axes of a central galaxies than the 
minor axes (jBr aincrd 20Q3). This, together with the sug- 
gestion that satellites s how a preference for radial align- 
ment with their hosts (|Agustsson &: Br ainerdl [20061 ). in- 
dicates that the measured flattening from galaxy-galaxy 
lensing may be underestimating the true halo flattening. 
The signal along the major axes will be suppressed by the 
intrinsic alignment of the satellites, thus our measure- 
ments should be considered a lower limit on halo elliptic- 
ity. In order to quantify this effect it would be interesting 
to measure the preferred alignment of satellite galaxies, 
and its dependence on their position angle with respect 
to their host galaxy for a large sample of hosts and satel- 
lites. Galaxy-galaxy lensing measurements using redshift 
information can eliminate this potential contaminant as 
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Fig. 11. — Ratio of mean shear experienced by sources closest 
to minor axis of a lens to those experienced by sources closest to 
the major axes as in Figure [9] only including just the lenses with 
axis ratio (b/a) between 0.5 and 0.8. This selection criteria will 
result in a lens sample with is dominated by early type galaxies. 
The weighted average shear ratio is 0.61±0.14. 

well as identify a sample of isolated foreground galaxies 
to eliminate the influence of the group or cluster envi- 
ronment. 

5. SUMMARY 

We have used early %' data from the CFHTLS to de- 
tect a significant galaxy-galaxy lensing signal, and to con- 
strain the velocity dispersion of galaxy halos at a redshift 
of 0.3. The measured velocity dispersion for an L* galaxy 
is 137±11 km s _1 (corresponding to circular velocity 
V* = 194±15km s" 1 ). This result is c onsistent with previ- 
ous galaxy-galaxy lensing e stimates (Klcinhci nrich et alj 
[2fM[Hoek"stra et all 12003) . 

In addition, we were able to estimate the total mass 
of an L* galaxy at redshift 0.3 to be 2.2±0.4 xlO 12 
ft _1 M0, which corresponds to a virial mass, M200 of 
1.1±0.2 xlO 12 /i _1 M Q . This is in good agr eement 
with the results from RCS jHoeks tra et all 120051). SPSS 
(iMandelbaum et al.ll2006al ). and GEMS (|Hevmans et all 
20063). 

We also saw some evidence for non-spherical dark mat- 
ter halos, though a definitive answer awaits more data. 
With multi-color data we will be able to repeat the 
anisotropic weak lensing measurements for sub-samples 
of lenses selected by color, morphology and redshift, This 
will also allow a detailed comparison to the low redshift 
results obtained by Mandclbaum et al. (2006a). 

Galaxy-galaxy lensing is a unique and powerful tool 
for studying the dark matter halos of galaxies to large 
projected distances. This technique has now been suc- 
cessfully applied to many different data sets, but with 
planned and ongoing major surveys the results will only 
improve. Any complete theory of galaxy formation must 
be able to to properly map both the baryonic and dark 
matter components of galaxies, and understanding which 
galaxies live in which halos is a major piece of this puzzle. 
Many observations can help in understanding the bary- 
onic component of galaxies, while simulations are best at 
explaining the dark matter. Gravitational lensing has an 
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important role to play in uniting simulations and other 
observations by connecting galaxies to their halos. 

There is clearly much more that can and will be done 
with the entire CFHTLS data set. In particular, having 
photometric redshifts will aid the analysis immensely. At 
present, we have lenses and sources over a wide range of 
redshifts and are estimating shear in angular bins which 
mix together many physical scales. Furthermore, pho- 
tometric redshifts will permit a clear division between 
lens and source populations, rather than a statistical 
one based upon observed magnitudes as used here. This 
should help alleviate the concerns of intrinsic alignments. 
Redshift information also permits the division of lenses 
into different redshift bins, thus allowing the study of 
halo evolution. 

A further use of the photometric redshifts could be 
to create catalogues of isolated lens galaxies. If the lens 
sample includes both central galaxies and their satellites, 
and the satellites are aligned with their host, then the 
galaxy-galaxy lensing signal could be biased. Also, the 
lensing signal around satellites will be influenced both 
by the dark matter halo of the satellite and the cen- 
tral galaxy. This issue can be addressed by measur- 
ing the lensing around only around the isolated "cen- 
tral" galaxies or by attempting to measure the signal 
around the central and satellite galaxie s separately (e.g., 
iHoekstra et"alll2005t lYang et al.ll2006ft . It is also pos- 
sible to mo del the satellite cont ribution to the lensing 
signal (e.g.. iGuzik fc SeliakJl27)0l . 



We also plan to measure the morphologies of all of 
the lens galaxies so that we can divide the sample into 
early and late-types and try to detect differences in their 
halo masses, shapes and M/L ratios. Furthermore, with 
the complete data set we should have sufficient statisti- 
cal power to be able to distinguish between dark mat- 
ter profiles. For example, we should be able to tell 
whether the tangenti al shear is best fit by an NFW profile 
(iNavarro et al.lll996l). an isothermal sphere or a Moore 
profile (jMoore et al.lfl999l ). 

We would like to thank Joel Primack and Tereasa 
Braincrd for helpful discussions. The authors would also 
like to thank the anonymous referee for a careful read- 
ing of the manuscript and many helpful suggestions that 
improved the paper. 
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Hawaii. This work is based in part on data products 
produced at TERAPIX and the Canadian Astronomy 
Data Centre as part of the Canada-France-Hawaii Tele- 
scope Legacy Survey, a collaborative project of NRC and 
CNRS. 
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